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We study the virial relations for ultracold trapped two component Fermi gases^ in the case of 
short finite range interactions. Numerical verifications for such relations are reported through the 
BCS-BEC crossover. As an intermediate step, it is necessary to evaluate the partial derivatives of 
the many body energy with respect to the inverse of the scattering length and with respect to the 
interaction range. They are found to have extreme values at the unitary limit. The virial results 
are used to check the quality of the variational wave function involved in the calculations. 

In the absence of interaction, the virial theorem relates the energy per particle of a confined atomic gas with the 
trapping potential. If that potential is harmonic, the theorem states that the total energy per particle is twice the 
mean trapping potential energy 

E = 2E tr . (1) 

For strongly interacting two component Fermi gases, confined by a harmonic trap in the unitary limit, this relation 
was also shown to be valid experimentally and theoretically^. The first derivation of that theorem considered zero- 
range interactions and made use of the local density approximation. Further insight on the fundamental basis of this 
relation revealed several remarking features of the unitary gas such as its scaling properties^ or a mapping, using 
group theory, between the trapped and the free space problem^. Recently, the Hellmann-Feynman theorem was used 
to prove Eq. (Q]) at the unitary limits, and to generalize the virial relations for finite scattering lengths^. In fact, 
general confinement potentials and finite range interactions can directly be taken into account using a general virial 
theorem which can be stated as follows^: 

Consider a Hamiltonian for a system of N particles with arbitrary statistics: 

H = H' + U(f 1 ,...,f N ), (2) 

where H' and its domain depend on p parameters with length dimensions l\, ...,i p , on K and the mass of the particles. 
U(fi, ...,r/v) denotes a regular arbitrary function that allows the domains of H and H' to coincide, fi is the position 
vector for the i-th particle. Then, 

B-(u + \±r,w)-\±t™. (3) 

x 1=1 ' q=l H 

with E the total energy. 

For N particles confined by a harmonic trap, U = m ( UJ x\ x i\ 2 + UJ 't\ViV + w zl z i| 2 )/2 and Eq. ([3]) becomes: 



E = 2E t . 
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where Et r = (U) is the trapping potential energy. 

In the present article, we study 2N fermionic atoms in two equally populated hyperfine states (TV = iVj = N± = 165) 
confined by an isotropic three-dimensional harmonic trap of frequency to, and interacting through an attractive finite 
range potential V = — |Vb|e _r / rv . This potential is characterized by two parameters, its strength Vb and its range 
r v . When the kinetic energy of the atoms is low enough, the scattering length is a proper parameter to describe the 
interacting system. For a given number of s-wave bound states and a given r v , there is a one-to-one relationship 
between the strength of the potential Vb and the scattering length a. We consider the case where at most one bound 
state is admitted by the potential and find the ground state of the many body Schrodinger equation approximately, via 
a variational Monte Carlo calculation, for several scattering lengths a and short potential ranges ?* v < < Th = y/h/mui. 
We then study the behavior of the total, internal and trapping energy as a function of both length parameters a and 
r v to verify Eq. Q . The explicit expression of the Hamiltonian is 



N ^2 _,_ ^2 



= X, + W (4 + rl) + E Vui, , (5) 
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and the corresponding virial relation becomes 

r v dE 
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r v — constant 



a— constant 2a 8(1/ a) 



r v — co7istant 
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In the BEC side of the crossover the total energy E can become extremely large compared to the total energy E in 
the BCS side due to the contribution of the binding energy of the formed molecules. This fact increases the numerical 
errors in the evaluation of the derivatives in Eq.®. In order to isolate this two-body effect from many-body effects, 
we have found convenient to take into account the behavior of the free space binding energy as follows. The two body 
problem, 



El 
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V(r) 



<p{r) = £<p{r), n = m/2, 



(8) 



is analytically solvable for s-states, so that the scattering length is explicitly given by 

a = r v 77(C) 



-2r v 



2 MO 



log(C/2) - C 



(9) 



with ( = (2r v y/\Vo\m/h), C = 0.577215664901... is the Euler constant and J v and N u represent the Bessel function 
of the first and second kind of order v 1 respectively. This problem has the following bound states 



(p{r) = JVJ 



>\m/h 



(y), 



(10) 



where N is a normalization factor and y = £e r / 2r ^ . The boundary condition at the origin implies J Xs (Q = 0, so the 



(r ) 

corresponding cigcncncrgies v fulfill the equation 

2r v \j \ei r ^\m/h = x s 



That is, x s is determined by £ and 
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We shall work with zq < ( < Z\ with zq and z\ the first two zeros of the Bessel function Jo. Under these conditions, 
just one bound state is admitted for each positive scattering length a. Given a and r v and using Eq. ([9]), we can write 



x Q = x a (0=x Q (r 1 - 1 (^))=w(^-). 



(13) 



As a consequence, the ground state binding energy £q of the two interacting particle system in otherwise free space 
satisfies the equation 
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Thus, if we define 
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if a < 0, 
if a > 0, 



r v — constant 
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and 
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the virial relation, Eq. ([7]), reads 
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This expression is easier to verify numerically than Eq. I|7]). Notice that, from a dimensional analysis, equations similar 
to Eq. (fl~3|) can be expected to be valid for other forms of the potential. 

Approximate ground state eigcnfunctions for the Hamiltonian Eq. were obtained variationally. The trial wave 
functions used have the Eagles-Lcggett form 



*x EL = A ] A l [<X1 T , li)0(2 T , 2{]..4(N U Ni)\ 



(18) 



through the BCS-BEC crossover regime. In this equation, A denotes the antisymmetrizing operator to be applied to 
all fcrmions of each species and 



(r iT ,r ji) = (p(n,j)e 



where ip(rij) is the s-wave ground state solution of the trapped interacting two body problem 
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cp(r) = sip(r). 



(19) 
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The variational parameter Xel modulates the optimal shape of the atomic cloud. The evaluation of the mean value of 
the many body Hamiltonian, Eq.([5]), for the Eagles-Lcggett trial wave function was done using Monte Carlo techniques 
that take advantage of the structure of the function^. 

For weak interactions, that is for negative scattering lengths shorter than the mean separation between interacting 
atoms, lower variational energies are obtained using the length scaled ground state solution of the noninteracting 
problem (which is a product of Slater determinants) multiplied by a Jastrow correlation function 



i=l,N 



i=l,N 



ill 



(21) 
(22) 



The scaling factors (3 and Xj were taken as variational parameters. For the many body ground state calculation, the 
inputs of the Slater determinants are the single-particle eigenstates of the noninteracting trapped system <^°(r), and 
the set of quantum numbers {n} are chosen to give the lowest energy compatible with Pauli exclusion principle. The 
Fermi energy tp derived by this procedure can be used to define an effective Fermi wave number Uf = v^mep /h, 
whose inverse is a natural unit for measuring the scattering length in the many body problem. 

The potential ranges used to perform the calculations were in the interval ry/rh o e[0. 002, 0.015]. For each ry and 
several scattering lengths through the crossover, upper bounds of the energy E were obtained by optimizing the 
variational parameters Xel or Aj and (3 according to the trial wave function. The optimization of the numerical 
subroutines allowed to explore higher statistics with respect to previous calculations^ and yield the evaluation of 
the energy with improved accuracy. In Figures 1 and 2 we illustrate the obtained energies, the vertical size of the 
plotted points are comparable or higher than the numerical errors. Notice that the resulting curves for the energy 
dependence on the scattering length for a given range and on the range for a given scattering length show a soft 
structure that allows a numerical interpolation or even an analytical local fitting curve. These interpolations were 
used to numerically compute the derivatives necessary to verify the virial relations. The results are illustrated in 
Figures 3 and 4. A remarkable feature found for all the short interaction ranges explored is that both derivatives, 
8E/8(l/kpa) and 8E/8ry, get an extreme value at unitarity. 

In Figure 5, we show a comparison between the trapping energy curve as predicted by the virial relation, Eq. (|17[) and 
specific values of that energy evaluated directly from the variational functions for a potential range ry = 0.0025r/,. o . 
Notice that there is good overlap between variational and virial results all over the crossover. The width of the 
continuous curve corresponds to the numerical errors in its derivation. The main source of error for the trapping 
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FIG. 1: Many body energy E, Eq. (|15[) . as a function of the inverse of the scattering length for a potential range rv = 0.0025rh o 
The energy units correspond to Eifg the energy of the non interacting atomic cloud. The range of the potential is measured 
in the length scale established by the trapping potential r^o = y^h/mui; kp denotes the Fermi wave number. 



energies evaluated directly from the variational wave functions results from the non uniqueness of the variational 
parameters Xel and /3 that yield similar variational energies. In fact, the agreement between the trapping energy 
evaluated directly and using the virial relation can be used as an additional criteria to select those parameters. As 
a reference, in Fig. 5 the many body energy curve E{l/kFd) is also shown. The crossing of virial (mcuR ) and total 
energy E curves does not occur at unitarity as a finite range effect. 

Summarizing, we have studied the virial relations expected for a balanced mixture of two species fermionic trapped 
atoms interacting through a two parameter (intensity and potential range) attractive potential. We have applied those 
relations to the study of the ground state solution of this problem when it is approximately obtained using simple form 
variational trial wave functions. In this way, we are able to quantify the quality of our wave functions and, even more 
important, we can compare the virial relations of short finite range interaction versus contact interactions. Along 
the calculations, the partial derivatives of the energy E as a function of the scattering length and of the potential 
range were numerically evaluated. It was found that in all cases considered those derivatives get extreme values at 
unitarity. As a consequence, an accurate determination of the coefficient E/Eifg at unitarity using finite range 
potentials require an accurate extrapolation procedure. In a similar way, an accurate experimental determination of 
such coefficient does require an extremely precise realization of the l/kpa — > limit. 
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FIG. 2: Many body energy £ as a function of the potential range rv at unitarity, 1/fcpa = 0. The energy units correspond to 
Eifg, the energy of the non interacting atomic cloud. The range of the potential is measured in the length scale established 



by the trapping potential = yjhjmuj. 
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FIG. 3: Partial derivative of the many body energy E, Eq. (|15p . with respect to the inverse of the scattering length as a function 
of the inverse of the scattering length for a potential range rv = 0.0025rh o . The energy units correspond to Eifg, the energy 
of the non interacting atomic cloud; &f denotes the Fermi wave number. 
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FIG. 4: Partial derivative of the many body energy E with respect to the potential range rv as a function of the inverse of the 
scattering length for a potential range ry = 0.0025rho- The energy units correspond to Eifg, the energy of the non interacting 
atomic cloud; &f denotes the Fermi wave number. 
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FIG. 5: (Color online) Softened curves for the many body energy E and the mean value of twice the trapping energy < mui R > 
as predicted by the virial relation, Eq. (|17[l . The dots represent the value of the trapping energy obtained directly from the 
variational wave functions. The energy units correspond to Eifg, the energy of the non interacting atomic cloud; kF denotes 
the Fermi wave number. 



